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The minor axes of, and the normals to, the polarization ellipses that surround singular lines 
of linear polarization in three dimensional optical ellipse fields are shown to be organized 
into Mobius strips and into structures we call "rippled rings" (r-rings). The Mobius 
strips have two full twists, and can be either right- or left-handed. The major axes of 
the surrounding ellipses generate cone- like structures. Three orthogonal projections that 
give rise to 15 indices are used to characterize the different structures These indices, if 
independent, could generate 839,808 geometrically and topologically distinct lines; selection 
rules are presented that reduce the number of lines to 8,248, some 5,562 of which have 
been observed in a computer simulation. Statistical probabilities axe presented for the most 
important index combinations in random fields. It is argued that it is presently feasible to 
perform experimental measurements of the Mobius strips, r-rings, and cones described here 
theoretically. 
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I. 



INTRODUCTION 



We describe here the Mobius strips, and other structures, generated by the axes of the polariza- 
tion elhpses that surround singular lines of linear polarization (L lines) in three dimensional (3D) 
optical ellipse fields. These strips have four half- twists; they complete the trinity of two-half- twist 
Mobius strips that surround ordinary ellipses [1,2], and the one-half- and three-half-twist strips 
that surround singular lines of circular polarization (C lines) [3] . 

On an L line the polarization ellipse collapses to a line, the major axis of the ellipse, and the 
minor axis and ellipse normal become undefined (singular). All three axes — major, minor, and 
ellipse normal — remain well defined for the ellipses that surround the L line. In a plane S pierced 
by an L line a point of linear polarization, an L point, appears. The projections onto S of the 
minor axes and ellipse normals of the surrounding ellipses rotate about the point with winding 
number (net rotation angle divided by 2tt) I = ±1 [4 — 13]. Typical structures that surround L 
lines are shown in Fig. [T} 



FIG. 1: Typical structures surrounding L lines in a (here computed) random 3D ellipse field. Shown in (a) 
— (d') (in (e)) are the structures generated by the minor (the major) axes of the ellipses whose centers are 
located in a plane So on a small circle that surrounds the line, (a) 3D view of a Mobius strip that surrounds 
an L line shown floating above its projection onto So. (a') Enlarged view of the projections of the ellipse 
axes in (a) onto Sq. Along a counterclockwise path (1 — 12) the axis projections rotate through 360° in the 
retrograde, clockwise direction, and / = —1. Other Mobius strips have / = +1. (b) The strip in (a) viewed 
from above. For clarity, here and in (b'), (c), (d), and (d') only half of each axis is shown. Here and in (d) 
ellipse centers are shown by gray circles and their minor axes by short straight lines; axis endpoints that lie 
above (below) the plane of the circle of ellipse centers are shown by filled white (black) circles. As can be 
seen, ellipse endpoints and centers form a pair of interlocking rings with four crossings. (b') The closed 
Mobius strip in (a),(b) opened to better display the double twist structure. Here the ellipse centers form 
a straight line (the thick line) around which the axes of the surrounding ellipses (thin lines) rotate. The 
small white (black) circles mark the arbitrary first (last) ellipse, (c) 3D view of a second type of structure 
that surrounds L lines, shown fioating above its projection onto So. This structure, which we call a rippled 
ring (r-ring), is not a Mobius strip, (c') Enlarged view of the axis projections onto So of the r-ring in (c); 
the projections rotate through 360° in the same direction as the path, and / — +1. Other r-rings have 
/ = —1. (d) The r-ring in (c) viewed from above; the endpoints and centers form noninterlocking rings, (d') 
The closed r-ring in (c),(d) opened to display the ripple of the axes, which oscillate through two complete 
sinusoidal cycles. In a random field approximately 1/3 of all L lines are surrounded by Mobius strips, 2/3 
by rippled rings. The normals to the surrounding ellipses (not shown) also generate Mobius strips and 
r-rings with / = ±1. (e) The major axes of the ellipses in So generate cone-like structures, a-cones, an 
example of which is shown floating above its projection onto Sq; here / = -|-1, other a-cones have / = — 1. 
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Different types of L lines are characterized by the geometrical and topological indices of the 
structures that surround the lines. Topological index I = ±1 generates two distinct types of L lines 
[4, 7, 11 — 13]; the Mobius strips and other structures introduced here greatly expand the number 
of these lines. In later sections we use 15 different indices to characterize these structures; these 
15 indices, if independent, could generate 839, 808 different L lines. This number is drastically 
reduced to 8, 248 by the selection rules described later; in a database containing 10^ independent 
realizations of a simulated random ellipse field we find some 5,562 different L lines. The true 
number of different lines is likely to be significantly greater, because as discussed later, due to the 
existence of highly improbable configurations the yield of different lines grows so slowly with the 
number of realizations that even 10^ independent realizations may be insufficient to generate all 
allowed possibilities. 

The plan of this report is as follows. In Section II we describe our computer simulations and 
the numerical and analytical tools we use in later sections to study L lines, their Mobius strips, and 
related structures. In Section III we describe in detail the Mobius strips, r-rings, and a-cones (Fig. 
[T| that surround L lines, and introduce 14 new indices to characterize these new structures. In 
Section IV we discuss the selection rules that constrain the number of possible index combinations, 
and present statistical probabilities for allowed combinations. We summarize our findings in the 
concluding Section V. Throughout, as in [1,3], our approach is descriptive, a more mathematical 
treatment similar to [2] will be presented elsewhere. 

II. METHODS 

We briefly review here the methods we use to study L points on L lines. To a large extent 
these methods are the same as those used in [3] to study C points on C lines, and the reader may 
find it helpful to consult this reference for more details. Many of the concepts used here and in 
[3] have their origin in [1,2], and the reader may also find it helpful to consult these references. 

The principal axis system for the general polarization ellipse is the (here right-handed) orthogo- 
nal three- frame a, /3, 7, where a and /3 are unit vectors directed along the major and minor ellipse 
axes, and 7 is a unit vector directed along the ellipse normal; the positive end of 7 can be uniquely 
defined by a right-hand rule based on the rotation of the electric field vector E as it traces out the 
ellipse over an optical cycle. 

As noted above, at an L point on an L line there is only one principal axis, a, that is well 
defined. Normal to this axis is a plane, the principal plane here labelled Sq. In what follows, we 
consider the 3D structure of the ellipses whose centers lie in Sq on a small circle ctq that surrounds 
the L point. We take the z-axis of the fixed, orthogonal coordinate system to be along a, and the 
corresponding xy-axes to lie in Sq, the orientation of these axes in this plane being arbitrary. 

As the plane of observation E is tilted relative to Sq the 3D arrangement of the ellipses whose 
centers lie in S changes. For very small tilt angles these changes affect only the geometries of the 
structures, but not their topologies, or their statistical properties. However, as the rotation angle 
increases past some small, critical value that differs for each L point, both the topology and the 
statistics change importantly. A similar phenomenon occurs for points of circular polarization, C 
points, on C lines [3]. All results presented in Fig. [T] (and unless stated otherwise in all other 
figures) are for the case S = Sq. The complex set of transformations that occur when S is rotated 
away from So will be reported on separately. 

In general, a, and therefore Sq, make arbitrary angles with the L line itself, so that in moving 
along the line from one L point to another the coordinate system must be rotated in order to 
remain in the principal plane Sq of each L point. 
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a, P, and 7 can be calculated from the (here complex) optical field E using either 

a = Re(EVE-E), (la) 
/3 = Im(EVE-E), (lb) 
7 = Im(E*xE), (Ic) 

which is due to Berry [11, 13], or from the eigenvalues, Aj, and eigenvectors, Ui, i = 1,2,3, of the 
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3x3 real coherency matrix [14] Mij = Iie{E*Ej), i,j = x,y,z. 7 as defined in Eq. (Ic) 
the area of the polarization ellipse and goes to zero at an L point. At an L point the polarization 
is linear and E can be made be pure real so that also /3 vanishes at the point. In what follows, a, 
f3, and 7 for the surrounding ellipses are, without change in notation, always normalized to unit 
length, i.e. a =^ a/ \a\, etc. 

We study two computed 3D ellipse fields. The first is composed of a large number of linearly 
polarized plane waves with random propagation and polarization directions, and random phases 
[1]. This field is an exact solution of Maxwell's equations and serves as an important check on the 
structures found using the more convenient linear expansion described below. L lines in this field 
were traced out using the L point discriminant Dl = 02 obtained from the characteristic equation 
A2 + aiA + a2 = of A% [15]. 

In the immediate vicinity of an L point the field describing the ellipses in So can be expanded 

as 



Ey — {Pyx -l- iQj/a;) ^ ~^ (^j/j/ ~^ '^Qyy) Vi 
Ez = a + {Pzx + iQzx) X + {Pzy + iQzy) y, 



(2a) 
(2b) 
(2c) 



where the direction of polarization of the point is along the z-axis and So is the xy-plane. 

In many cases simpler expansions suffice: the Mobius strip in Fig. [T^,b is closely approximated 
by the field model Ex = ix + {1 — i)y,Ey = —iy,Ez = —1, the r-ring in Fig. by Ex = 

— {1 + i) X, Ey = ix + {1 — i)y, Ez = 1, and the a-cone in Fig. [l^ by E'a; = —x, Ey = x — y, Ez = 1. 

The statistics in Sq of a and the P, Q in Eq. ^ were obtained using the numerical procedure 
described in [3] for C points, modified for L points by the fact that whereas E • E = for a C point, 
E*xE = for an L point [11 — 13]. The probability density functions (PDFs) in Sq of a, and of 
P, Q are shown in Fig. [2] — here all P, Q have the same PDF. 
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FIG. 2: PDFs in So of the parameters a,P, and Q, in Eq. (a) PDF of a. (b) PDF of P,Q. 

numerical data (histogram) closely matches a Gaussian with unit variance (curve). 
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A single projection, the Sq projection, cannot, of course, fully describe the complex 3D arrange- 
ment of the ellipses surrounding an L line, and in addition to projecting the ellipses on ctq onto Sq 
we study two additional projections — tq and ttq. These projections, which are the same as those 
used in [1 — 3], are reviewed in Fig. [3| Thus, in characterizing the 3D arrangement of the ellipses 
in the vicinity of an L line we use three orthogonal projections, Sq, tq, and ttq — the minimum 
number required for a complex object. 




FIG. 3: Projections tq and ttq. The L point is at the origin of the fixed xyz-coordinate system, (a) Rotating 
plane tq, the XZ-plane {ttq, the FZ-plane), is radial (tangential) to ctq and travels on ctq, rotating about 
the fixed z-axis; its position on ctq is measured by the angle Xi and by arc length s — tqX- An axis, a, 
f3, or 7, of the ellipse whose center (small gray circle) is located at the origin of XYZ is projected onto tq 
(ttq). The endpoint of the axis is shown by a small white circle, the axis itself by a thick line connecting the 
center and the endpoint. The endpoint projection in tq (in ttq) is shown by a small black circle. The vector 
T (P) in To (in ttq) connects the ellipse center with the endpoint projection. As tq (ttq) moves along ao 
different ellipses are projected onto tq (ttq), their endpoints in tq (ttq) trace out a curve, the tq (ttq) endpoint 
curve, and vector T (P) rotates. As discussed in Section III, the winding numbers of rotation of T (P) and 
of the tangent to the endpoint curve in tq (in ttq), give rise to indices Ti and dxi (tt^ and d'Ki),i = a,/3, 7, 
that characterize the 3D arrangement of the ellipses on ctq- (b) Special projections. When the endpoint 
projection onto tq (onto ttq) of an ellipse axis lies on the Z-axis, as shown here for an axis of ellipse 1 (ellipse 
2), the axis projection onto Sq is radial to (Tq, i.e. along the X-axis (is tangential to ctqj i-e. along the 
y-axis). These two important special cases are elaborated on in Section III. 

Issues regarding scaling and other technical aspects of the L point graphics presented here are 
the same as those for the C point graphics presented in [3] ; the reader will find a full discussion of 
these issues in [3] . 

Unless stated otherwise, throughout we use a single value for the radius tq = 10~^ of the 
surrounding circle ao. This value is sufficiently small that except for a scale factor the wavefield 
structure is invariant under further reduction of the radius. The fact that every small value for rg 
yields the same structure implies that the L point is surrounded by nested Mobius strips, r-rings, 
a-cones, etc. The full 3D arrangements of these deeply nested structures are, of course, quite 
impossible to visualize; accordingly, below we dissect out a single, typical surrounding circle, the 
Ctq circle with ro = 10~^, and proceed to study its properties in detail. 
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III. INDICES OF MOBIUS STRIPS, R-RINGS, AND a-CONES 

Here we discuss the 12 winding numbers that describe the Mobius strips, r-rings, and a-cones 
that surround L hues, as well as the line classification [16] and its relationship to these indices. 



A. Indices of the Projection onto Eq^ ^/3,7 /, 



As a result of the studies of Nye and coworkers [4 — 10] , it has long been known that the winding 
number of the projection of axis /3 onto Sq, Ip = ±1 (Figs. [l^',c'), equals the corresponding 
winding number of axis 7. The reason for the equivalence of these two indices is that as the radius 
of the surrounding circle shrinks towards zero, by continuity the major axes of the ellipses on this 
circle all approach parallelism with the direction of polarization of the L point. The projections 
of j3 and 7 onto Sq are therefore substantially orthogonal, and as one rotates so does the other. 
Accordingly, in what follows we use the symbol Ip^^ to denote the common winding number of 
these two axes. Although perhaps less obvious, the equivalence I 13 = holds also for arbitrary S 
[13]. 

The winding number of the projection of axis a onto Sq does not appear to have been 
discussed previously. We find that just like for 1^ ,^, in all cases = ±1 (Fig. [l^). 

la and Ifs^^ are independent, and we find all four possible sign combinations in our simulations. 

It is not difficult to derive analytical expressions for the above indices. Writing V for axis cx, 
j3, or 7, we have from Eqs. ([T]) and ^ for the j;y-components of V for sufficiently small x, y, 



— PxxX ~\~ PxyVi 



(3a) 
(3b) 



where for 



axis OL 
axis 13 
axis 7 



F = P (i.e. F,, 

^ XX - 



Pxxi etc.). 



lOjQyX , Fxy 



IdQyy^ Fyx 



(4a) 
(4b) 
(4c) 



Berry [13] has given the formula / = sigii{FxxFyy 
form in Eq. (|3|. Inserting Eq. (|4| into Berry's formula, we have 



FxyFyx) for a general vector field of the 



la — ^'^S^ (.-^xxPyy PxyPyx) > 



'/3 



-^/3,7 — sign (^QxxQyy QxyQyx^ • 



(5a) 
(5b) 



We find this result to be in full agreement with the indices obtained numerically. 
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B. Indices of the Projections onto tq and ttq 



1. Projections onto tq 



a. Indices Ti and dri, i = /3, 7 



Mobius strips Both axes (3 and 7 of the surrounding eUipses can generate Mobius strips 
that are either right- or left-handed screws with two full turns. 

Index Ti measures the number of turns and their handedness (sign), and is Ti = +2 (Ti = —2) 
for a left-handed (right-handed) Mobius strip. As can be seen, the Mobius strip in Fig. [T^— b' is 
a left handed screw with = +2. 

Although not an obvious geometrical or topological necessity, we find in Sq that if rg = +2 
(rg = —2) then 7^ —2 (r^ 7^ +2), and vice versa. If both axes /3 and 7 generate Mobius strips 
(which need not, and does not, always occur), then the above rule implies that both strips must 
have the same handedness. 

Index dT measures the rotation of the tangent to the endpoint curve. For Mobius strips, and 
only for Mobius strips, we find that in all cases dTi = Ti, i.e. when Ti = +2 (Ti = —2) dTi = +2 
{dTi = —2). Because Tp = r-y when both axes generate Mobius strips, dTi = Ti implies that for 
such paired strips drg = dr^. 

Indices Xj and dTi for L point Mobius strips are illustrated in Fig. |4j 



Z 
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FIG. 4: Mobius strip endpoint curves in tq. Here and throughout, small, connected black circles show the 
curve generated by the projections of axis endpoints onto tq as tq rotates on ctq through one complete circuit 
around the L point; the starting point of the circuit is shown by a small white circle; the sequence in which 
the endpoint curve is drawn is indicated by a small arrow. (a),(b) tq endpoint curve of the Mobius strip in 
Figs. ^-h\ (a) 3D view of the endpoint curve. The first and last points are shown connected because they 
are the same point. The central straight line is formed by the centers of the ellipses on ito- As can be seen, 
the endpoints generate a two-turn left-handed helix, (b) The endpoint curve in (a) viewed normal to tq. The 
"x2" label indicates that the second turn of the helix in (a) lies behind the first turn visible here. Vector 
T has its origin at the center of the ellipse on the surrounding circle co (the origin oi XZ) and its head at 
the endpoint projection. The winding number r^, where for the Mobius strip shown here i = (3, is obtained 
from the rotation of T by measuring the rotation angle i^{s) — arctan (T^ (s) ,Tx (s)), and unwrapping 
(unfolding) the result to obtain Ti — [il> {2tt) — V (0)] / (27r). The winding number dr^ of the tangent vector 
to the endpoint curve, dT, is similarly calculated from winding angle aictan (dTz /ds,dTx /ds). Here the 
net total rotation of T and of dT is +47r, and Ti — dTi — +2. Most endpoint curves do not have a 
symmetrical, near circular shape, as does the one shown here, which was chosen for simplicity and clarity 
of presentation. Although here = dr^, in general these two indices can be different. 
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An important property of L point Mobius strips is that the first and second turn of the two-turn 
hehx are identical, so that the projections onto tq of the two turns always overlap and only one 
turn is seen in such a projection, Fig. [4| We therefore include in such projections a "x2" label to 
alert the reader to the fact that there is a second turn underlying the first. 

r-rings The endpoint curves of r-rings form two-turn helices — possibly surprising because 
such a structure cannot be easily inferred from Figs. [TJ;— d'. For r-rings, Tj = 0, whereas dri = 
0, ±2. Typical examples are shown in Fig. |5] which also illustrates the "phase ratchet rules" that 
simplify calculation of for complicated endpoint curves [3] . For paired /3 — 7 r-rings drp and 
dr^ need no longer be equal, and we find in our simulations all nine possible combinations of these 
two indices. 




FIG. 5: r-ring endpoint curves in tq. (a) The opened r-ring in Fig. [Tji' viewed from the left to display 
the two-turn helix generated by the axis endpoints. (b) Endpoint curve of (a). Because the endpoints of 
an r-ring do not spiral around the line of ellipse centers, the origin of the vector T always lies outside the 
closed endpoint curve, and as a result Ti = 0. This can be seen by following the rotation of T or by using 
the phase ratchet rules. At turning points ti and t2 vector T reverses its direction of rotation, which is 
counterclockwise, i.e. positive, (clockwise, i.e. negative) on the right (left) half of the endpoint curve, and 
the net rotation of T is zero. In applying the phase ratchet rules [3], one first lists the signed crossing of the 
endpoint curve with the XZ-axes in the order that they occur, starting at the small white circle that marks 
the beginning of the circuit on (Tq and proceeding in the direction of the small arrow. For the endpoint 
curve of the Mobius strip in Fig. |4]d the crossing sequence is Z+X+Z+X^, where the sign of a crossing is 
positive (negative) if T rotates counterclockwise (clockwise) at the crossing. This is the canonical sequence 
for Ti = +1] multiplying by 2 to account for the hidden second turn yields the net index — +2. For 
the r-ring shown here the crossing sequence is Adjacent terms in the same axis necessarily have 

opposite signs and cancel, so for this sequence n = 0. The tangent vector, dT, however, is easily seen to 
rotate through 27r in the positive, clockwise direction; multiplying by 2, the net winding angle is 47r and 
dTi = +2. (c) An r-ring whose endpoint curve forms a figure of eight. Here the crossing sequence is X_|_X_ 
and once again Tj = 0, as it is for every r-ring. But here also dr^ = 0, as it is for every figure of eight. That 
the net rotation of dT is zero can be seen by following the tangent vectors marked with the cumulative 
rotation angle, (d) An r-ring with = and dr^ = —2. In a random field there are roughly equal numbers 
of each type of r-ring: = and dri = —2, 0, +2. 



9 



b. Index dT^ As illustrated in Fig. [T^, axes a of the surrounding ellipses are not organized 
into Mobius strips or r-rings, but rather into cone-like structures — a-cones. In Fig. [6] we show 
additional views of the cone in Fig. [T^. As can be seen, the endpoint curve of an a-cone is a two- 
turn helix, much like the helix of an r-ring, and like an r-ring, Tq = also for the cone; once again 
the reason is that the endpoint curve does not encircle the ellipse centers and vector T therefore 
oscillates clockwise and counterclockwise with a net winding angle of zero. Unlike r-rings, however, 
the endpoint curves of a-cones are never figures of eight, and for the cones dTa = ±2. 

For the a-cone in Fig. [6^ the axis endpoints form an ellipse that is smaller than the circle of 
ellipse centers (i.e. the surrounding circle). But for other cones, or for this cone viewed from below 
rather than from above, the endpoint ellipse is larger than the circle. In either case the two curves 
— circle and ellipse — never cross. There are, however, numerous examples of cones where the 
major axis of the endpoint ellipse is larger than, and the minor axis is smaller than, the diameter 
of the circle of ellipse centers, and in these cases the two curves — circle and ellipse — cross four 
times. 




FIG. 6: a-cones. Shown are three additional views of the cone in Fig. [Tp. (a) The cone viewed from 
above, (b) The cone in (a) opened and straightened to show the osciUations of the axis endpoints. (c) 
The To projection of the axis endpoints. Here the Z-axis is elongated relative to the X-axis to show the 
structure of the normally very flat two-turn endpoint helix that is characteristic of a-cones. 



2. Projections onto ttq 

a. Indices vTj and dni, i = (3,^ Indices vTj and di^i are the vro projection analogs of indices Tj 
and dTi of the tq projection. 

Index vTj measures the rotation of the endpoint curve in ttq, whereas index dvrj measures the 
rotation of the tangent to this curve. Fig. |4] is applicable to vrj and dvrj upon replacement of the 
X-axis by the y-axis, and vectors T and dT by vectors P and dP. ttj and diii can each take on 
the values 0, ±2. 

\Tp\ = 2 (|r-y| = 2) implies that axis j3 (axis 7) generates a Mobius strip, because for r-rings 
(and cones) Tj = 0. But the same is not true for vTj which equals ±2 if the endpoint curve in ttq 
encloses the ellipse centers, and is zero otherwise. In very nearly 1/4 of all r-rings |7rj| = 2. 

There exists a vTj selection rule that is the analog of the Tj selection rule noted above: when 
TT^ = +2 (tt^ = -2), / -2 (tt^ / +2). 

dvTj is ±2 for simple (non-self-intersecting curves), and is zero for figures of eight. When vTj = +2 

(vTi = -2), dT^i = +2 (dvTi = -2). 

Fig. [7] illustrates the case of an r-ring with Tp = drp = 0, and vr^ = dvr^ = —2. 
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FIG. 7: r-ring endpoint curves of axis 7 in tq and in ttq. (a) tq endpoint curve. Both and dr-y are zero, 
as is always the case for a figure of eight, (b) ttq endpoint curve. The curve encloses the origin (circle of 
ellipse centers), and both and dn^ equal —2. 

b. Index d'Ka Indices vr^ and dn^ are the ttq projection analogs of indices Tq, and dra of the 
To projection. Like the endpoint curve of axis a in tq, for which = and dTa = ±2, the 
endpoint curve in ttq generates a two-turn helix for which tTq = and dua = ±2. Fig. |6j; for the 
To projection is applicable also to the ttq projection after the X-axis is replaced by the y-axis. 

This completes our discussion of the 12 indices that characterize the three orthogonal projec- 
tions, So, To, and vro, of the Mobius strips, r-rings and a-cones that surround L lines. In the next 
section we return to the So projection and consider a characterization [16] that is well known for 
C lines [4, 6, 11 — 13], but that does not appear to have been discussed previously for L lines. 



C. Line Classification: Indices Ai,i = a,/3,7 
1. Computer Simulation 

The arrangement of the axis projections onto So can be characterized by two indices. The 
first, winding numbers Iq, and Ip = = Ip^^, measures the net rotation of the axis projections 
around the L point, and has been discussed in Section III. A and illustrated in Fig. [T] The second, 
the line classification index Ai,i = a, f3, 7, counts the number of streamlines formed by axis i that 
terminate (or originate) on the L point as straight lines (separatixes). For a sufficiently small 
surrounding circle do, as is used here, Aj counts the number of axis projections on o"o that point 
directly at the L point. We find that in all cases, for positive (negative) /j. A, = 0,4 (4). We 
note that although the line classification [16] has not been previously discussed explicitly for L 
points for any axis, these values are consistent with streamline maps for axes f3 and 7 presented 
in [4, 7, 11]; streamlines for axis a do not appear to have been considered previously. 

Unlike the case of winding number I which is always the same for both axes /3 and 7, A^ 
and A^ can differ when Ip^^ = +1. Examples of the line classification are illustrated in Fig. [8| 
An important, general property of these maps is that they posses a center of inversion. This 
follows from the forms of Eqs. ([T]) and ^ for sufficiently small x, y. But this is the only intrinsic 
symmetry of such maps, and therefore of the corresponding streamlines. 

As indicated in Fig. |3] an axis projection points towards the L point (circle center) only when 
the axis itself lies in the XZ-plane, so that the projection of its endpoint onto vro lies on the Z-axis. 
Also this property is illustrated in Fig. [8} Thus, A also equals the number of times the endpoint 
curve in ttq crosses the Z-axis, and vice versa. 
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FIG. 8: Line classification in a computed 3D random ellipse field, (ai) — (as) Axis projections onto Sq. 
Shown as short straight lines are the projections of the axes of the ellipses on the surrounding circle ctq- 
Projections that point directly at the circle center (L point), here A- lines, are marked by thick lines; the 
line classification index A counts the number of these lines. In our simulation A-lines are found numerically 
[3] as the solutions of 9 {x) = X + "■"'i with n — —4. ..4, as required. Here is the angle that the axis 
projection makes with the x-axis, and x measures the position of the ellipse center on a^, Fig. |3] The 
straight lines passing through the L point are calculated from Eq. As can be seen, in each case shown 
here, and in numerous other cases tested but not shown, there is complete agreement between theory and 
the simulated field, (ai) Axis a. Here index = —1, Aq. = 4; for other a-cones 1^ = +1, Aq = 0,4. 
{a.2) Axis (3. Here //3 = +1, A^ = 0; for other (3 Mobius strips or r-rings = — 1,+1, A^ = 4. (aa) 



Axis 7. Here 



-1, A^ = 4; for other 7 Mobius strips or r-rings 



-1, A. 



4, or 



+1, 



A^ — 0. (bi) — (bs) Axis projections onto Sq corresponding to the ellipse fields in (ai) — (aa), respectively. 
Shown are axis projections for ellipses whose centers lie on a square grid that surrounds the L point. These 
projections (short lines) are the tangents to the streamlines; visually, their overall pattern approximates the 
streamline pattern. The long straight lines are the separatixes calculated from Eq. ([9]), and as can be seen 
are in full accord with the patterns shown here, as well as with numerous other examples of these patterns 
that are not shown, (bi) This streamline pattern is similar to the contour pattern surrounding a saddle 
point. (b2) These streamlines spiral around the L point. ih^) Here the streamlines form right and left 
leaning parabolic-like curves "centered" on the right leaning separatix; the fields of right and left leaning 
"parabolas" are separated by the left leaning separatix. These three patterns are the only ones found in 
our simulation. We find that each pattern is tied to a particular combination of / and A: in all of the 
numerous cases tested for axes a., (3, and 7, when / = —1, A = 4, the pattern is qualitatively similar to that 
in (bi); when I = +1 and A = the pattern is qualitatively similar to that in (b2); and when / — +1 and 
A = 4 the pattern is qualitatively similar to that in (ba). (ci) — (C3) Endpoint curves in ttq corresponding 
to the axis projections in (ai) — (aa). As can be seen, and as is expected (Fig. |3|, A equals the number of 
Z-axis crossings of the endpoint curve in ttq. This has been verified for all 10^ realizations in our random 
field database. 
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2. Theory 

The above results can be obtained analytically from the properties of the streamlines. These 
streamlines, expressible as the family of curves y (x), are the solutions of 



dy^V^^ FyxX + Fyyy 
dx Vx FxxX + Fxyy 

where in writing the second form we have used Eq. Fxx, Fxy, etc., are given in Eq. Q. 

Eq. ([g]) is a classic exercise in the theory of ordinary differential equations, and the properties 
of its solutions are well known [see for example 17]. There are six types of solution, three of 
which are unstable under small perturbation, and are therefore nongeneric and do not appear in 
our simulations, and three which are stable and which do appear. For completeness, and in order 
to help the reader avoid a possible misinterpretation of previous illustrations of /3,7 streamlines 
[4, 7, 11], below we consider briefly all six types. 

The solutions of Eq. mn are characterized by both index / and the discriminant 



A — {Fxx — Fyy)^ + AFxyFyx- (7) 

When A < 0, A = 0, when A = 0, A = 2, and when A > 0, A = 4. As discussed below, the 
case A = 0, A = 2 corresponds to one of the unstable solutions and is not seen in our simulation; 
the other two possibilities are. 

The rule discussed above, that for Ij = —1, Aj = 4, i = a, /3, 7, can be easily proven by rewriting 
Eq. ^ as 

^ {Fxx ~l~ Fyy) 4X, 

FxxFyy FxyFyx- (8) 

Recalling from Eq. ^ that I = sign(X), it immediately follows that if / < 0, A > 0, and as a 
result, A = 4. 

The six solutions of Eq. ([g]) are illustrated in Fig. [9| where the three unstable (stable) solutions 
are labelled Ui, U2, and U3 (Si, S2, and S3). 

Ui corresponds to dy/dx = y/x, and requires Fxx — Fyy and Fxy — Fyx — 0. For these values 
of the parameters A = 0, and for this special case h. = 00 rather than A = 2. When the stringent 
conditions on F are no longer met, A 7^ 0, and Ui transforms into the stable configuration S2 (S3), 
if A < (A > 0), in both cases preserving index /, as expected. 

U2 corresponds to the condition Fyy = —Fxx- When this special condition is met the streamlines 
form conic sections: ellipses if A < 0; hyperbolas if A > 0. These hyperbolas are organized into 
the contour pattern of a saddle point. When the special condition on F is no longer met the 
ellipses open into the stable spiral S2, whereas the pattern of the hyperbolas is maintained as the 
stable solution Si, although the curves themselves are no longer hyperbolas. As expected, in both 
cases index I is conserved during the transformation. 

U3 corresponds to the special condition FxyFyx = —\{Fxx — Fyy)"^. When this special condition 
is met A = and A = 2. When the condition on the F' s is no longer met U3 transforms into the 
stable solution S2 (S3) if A < (A > 0); again, index / is conserved during the transformation. 

Worth noting is that when the special conditions on F are very nearly satisfied, as can and 
does happen in a random (or other) field, the pattern of streamlines can closely approximate the 
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unstable solutions, possibly resulting in the misleading impression that these solutions have been 
observed, and are therefore stable. 

For the three stable solutions Si, S2, and S3, indices / and A are the following: for Si, I = —1 
and A = 4; for S2, / = +1 and A = 0; for S3, / = +1 and A = 4. These index combinations 
uniquely determine the streamline pattern, and parallel in some sense the three index combinations 
for the three types of stable C points, "star" , "lemon" , and "monstar" [4,11 — 13, 16] , although the 
patterns themselves differ for the two cases. The C point combinations are: for the star, I = —1 
and A = 3; for the lemon, / = +1 and A = 1; for the monstar, / = +1 and A = 3. 

When A = 4, the straight-line streamlines passing through the origin (L point) are given by 



y± 



yy 



2R 



-X. 



(9) 



xy 



These lines are separatixes that divide the field into four quadrants in which the sign of the 
curvature of the streamlines changes across a separatix. 

We find that our simulations (Fig. 8) are in full agreement with all of the above, and that all 
three stable solutions Si, S2, and S3, are present for each axis a, /3, and 7, and for each type 
of structure — a-cones, r-rings, and Mobius strips. This reflects the fact that for each of these 
structures all three index combinations, / = — 1,A = 4 (Si), / = 1,A = (S2), and / = 1,A = 4 
(S3), appear. 







FIG. 9: Solutions of Eq. U1-U3 (S1-S3) are unstable (stable) solutions. The parameters used 

in plotting these solutions are for the unstable solutions: Ui, F^x — ^,Fxy — Fy^ — 0,Fyy = 1; U2, 
Fxx ^ 1/2 - S {S ^ ^ S = O.l), Fxy = l,Fyx = -l,Fyy = -1/2; U3, F,, = l,Fxy - l,Fyx = 
— 1/4, Fyy = 0. For the stable solutions the parameters are: Si, F^x — l,Fxy — ~l,Fyx — l,Fyy — —2; S2, 
Fxx = l,Fxy = -2,Fyx = S,Fyy = 4:] S3, Fxx = ^ , Fxy = 1/2, Fyx = l,Fyy = 2. THese parameters were 
chosen to yield plots that approximately match the streamlines in Fig. [8] Si approximately corresponds to 
the pattern in Fig. S2 to the pattern in Fig. ^2', S3 to the pattern in Fig. 



D. Index Summary 



We briefly summarize here the major properties of the 15 indices that characterize the Mobius 
strips, r-rings, and a-cones that surround L lines. 
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a. la, Aq,, Ap, A^. These five indices characterize the projections of axes a, /3, and 7 
onto the principal plane Sq. Index / measures the winding number of these projections, A counts 
the number of straight streamlines (separatixes) that terminate (or originate) on the L point. For 
all three axes I takes the values ±1 and A takes the values 0,4. If I = —1, A = 4, whereas for 
/ = +1, A = 0,4. Each of the three possible combinations of / and A defines a unique streamline 
pattern. Figs. [8] and [9} 

b. dTa, T/3, dTi3, T^, cIt^. These five indices characterize the projections of axes a, /3, and 
■~f onto the plane tq. Fig. |3j r measures the winding number of the axes about the circle of 
ellipse centers, ao- For Mobius strips r = ±2, where +2 (—2) implies a left-handed (right-handed) 
two-turn helix. Figs. [T]and|4| For r-rings (and a-cones) r = 0. 

dr measures the winding number of the tangent to the endpoint curve. If this is a simple closed 
curve dr = +2 {dr = -2), Fig. |4]and|6j If the endpoint curve generates a figure of eight, which 
can occur for axes /3 and 7, but not for axis a, dr = 0, Fig. [5] When r = ±2, dr = r. 

c. diTa, T^j3, dTTjS, TT^, diT^. These five indices characterize the projections of axes a, /3, and 7 
onto the plane ttq. Fig. [sj If the endpoint curve in ttq encloses the origin (circle of ellipse centers), 
which can occur for Mobius strips and r-rings, but not for a-cones, vr = ±2, otherwise vr = 0. dir 
measures the rotation of the tangent to the vro endpoint curve. The possible values of this index 
for the different axes are the same as those for dr. When vr = ±2, dn = n. Nz, the number of 
times the endpoint curve in vro crosses the Z-axis, is equal to A, the line classification index, so Nz 
could replace A; no index or combination of indices in ttq and/or tq, however, can replace /. 

For a given structure — Mobius strip, r-ring, or a-cone — indices r and dr of tq, and indices 
TT and dn of ttq, can be, and often are, different, and to a large extent these indices are also 
independent of the indices / and A that describe the projection onto Sq. There are, however, 
index combinations that do not occur; in the following section these missing combinations are 
summarized in the form of selection rules. 

IV. STATISTICS 

We discuss here the probabilities of various combinations of the 15 indices that characterize 
the Mobius strips, r-rings, and a-cones that surround L lines, starting with those combinations 
for which the probability is zero. These zero probability combinations are summarized below in 
the form of binary and ternary selection rules. Quaternary and higher-order rules may also exist. 
Discussed below are simple, systematic methods for identifying all true binary and ternary rules; 
in contrast, there do not appear to be equivalent methods for identifying true higher-order rules 
(i.e. rules that do not incorporate combinations of lower-order rules). 

A. Selection Rules 

The binary and ternary selection rules listed below were either derived from results discussed 
earlier, or deduced from the absence of various configurations in our 10^ entry database. These 
latter, empirical rules await confirmation from theory. 

Finding the empirical rules was facilitated by the use of 2D and 3D correlation (scatter) plots. 
Of the 105 2D plots of two index combinations, 17 had missing configurations that led to the binary 
selection rules listed below. The number of 3D, three-index plots was 455; of these 198 had missing 
configurations, most of which, however, involved combinations of missing two-index configurations 
already included in the binary rules. The ternary selection rules listed below were obtained from 
the 12 plots that contained true three-index missing configurations (some plots yielded more than 
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one rule) . Typical examples of 2D and 3D correlation plots with missing configurations are shown 
in Fig. [T0| 




FIG. 10: Index correlation plots. The relative areas of the circles approximate the relative probabilities of 
the different configurations, (a) 2D plot of the index combination Tfj — t^. Empty grid points correspond 
to missing configurations, (b) 3D plot of the index combination — Tp — t^. Missing configurations are 
marked by x's. The reader may enjoy identifying the selection rules obtained from these plots. 

The binary (ternary) selections rules acting alone reduce the number of configurations of all 
15 indices from 839, 808 to 16, 344 (98, 864); both sets of rules acting together yield 8, 248 allowed 
configurations, of these 5, 562 have been observed in our simulations. 

1. Binary rules 

These rules involve pairs of axes. They state that the following two-index configurations are 
forbidden: 



h 




1 & Aj = 0, i 


= a,/3,7- 


(Rule 2.1) 




= 2 


& Ai = 0, i = 


/3,7- 


(Rule 2.2) 


\Ti\ 
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(Rule 2.3) 
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(Rule 2.4) 


diTi 
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(Rule 2.5) 
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= 2 


& d-Ki / TTj, i 


= /3,7. 


(Rule 2.7) 




= 2 


& |r^| = 2 & 


sign (rg) / sign (r^) . 


(Rule 2.8) 




= 2 


& Ivr^l = 2 & 


sign (vT/j) / sign (vr^) . 


(Rule 2.9) 



The above rules either follow from already discussed, or imply new, geometrical connections 
between the structures formed by the different axes for a given set of wavefield parameters. Below 
we briefly discuss these connections. 

Rule 2.1. This rule is proven in Eq. ([s]). 

Rule 2.2. This rule follows from an already discussed geometrical connection: |7r| = 2 requires 
the endpoint curve in tto to cross the Z-axis four times. Section III. B. 2a and Fig. |4j at each such 
crossing the axis projection is radial. Fig. [3j and so A = 4. 

Rule 2.3. This rule follows from a an already discussed geometrical connection: |r| = 2 (i.e. 
a Mobius strip) requires the endpoint curve in tq to cross the Z-axis four times; at each such 
crossing the axis projection is tangential. Fig. [3} Because, as already noted in Section III.A, the 
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projections of axes /3 and 7 onto Sq are orthogonal, where the (3 (7) projection is tangential the 
7 (/3) projection is radial, and so = 4 (A^ = 4). 

Rule 2.4- This rule implies a new geometrical connection: if the endpoint curve in tq is a figure 
of eight, it must cross the Z-axis four times (see the explanation of Rule 2.3). 

Rule 2.5. This rule implies a new geometrical connection: if the endpoint curve of axis (3 (7) 
in TTo is a figure of eight, the endpoint curve of axis 7 {(3) in ttq must cross the Z-axis four times. 

Rules 2.6 and 2.7. These rules follows from already discussed geometrical connections: if 
|r| = 2 (|7r| = 2) the endpoint curve is elliptical in shape, and therefore dr = t {d-K = vr). Fig. [4j 

Rule 2.8. This rule implies a new geometrical connection: if both axis (3 and axis 7 generate 
Mobius strips, both strips must have the same handedness. 

Rule 2.9. This rule implies a new geometrical connection: if the endpoint curves in ttq of axis 
(3 and axis 7 are both elliptical in shape, both curves must have the same sign. 



2. Ternary Rules 



These rules involve three different axes. They state that the following three-index configurations 
are forbidden: 
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All the ternary rules involve connections between a-cones and Mobius strips or r-rings. These 
connections, all of which are new, appear difficult to fathom in easily understood, simple geometrical 
terms. For example. Rule 3.1 states that if the streamlines of axis a form a spiral. Fig. |9]-S2, and 
those of axis /3 (axis 7) form a saddle point. Fig. |9]-S1, then axis /3 (axis 7) itself must generate a 
Mobius strip. 



B. Probabilities 



The relative probabilities of all 15-index combinations found in our simulation are displayed 



in Fig. 11 in the form of a modified Zipf plot. In this plot the number of occurrences of a 



configuration, A^, is plotted vs. the rank of the configuration, R, where configurations are ranked 
in descending order of probability. As can be seen, for R > 1, N decreases step- wise to one with 
increasing rank in a quasi-continuous fashion. 
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But if N decreases quasi-continuously to one, then what meaning can be attached to the empir- 
ical selection rules for which N = 0? The answer is that there is a qualitative difference between 
the two-index (three-index) and the fifteen-index cases, because in the case of the two-index, binary 
(three-index, ternary), selection rules, the number of observed occurrences is either greater than 
10^ (10^) 



or zero. 



Shown in the inset of Fig. 11 is the number of unique index combinations, n, as a function 
of the number of realizations, r, where r is incremented sequentially in steps of 10^. As can be 
seen, n appears to asymptote to the power law n = 2, 244r'''^^^^. The parameters in this law 
were chosen such that for r = 10^, i.e. for our simulation, n = 5562, i.e. the number actually 
observed. The binary and ternary selection rules yield 8, 248 possible configurations. Assuming 
that higher-order rules do not significantly reduce this number, and that the power law remains 
unchanged for much larger r, some 4 x 10® realizations would be needed to capture all possibilities. 
More realistically, the yield of new configurations can be expected to fall off with increasing r, so 
that even 10^ realizations might not suffice to yield all allowed configurations. 




1000 2000 3000 4000 
Rank 



5000 6000 



FIG. 11: Modified Zipf plot. Shown is the number of configurations N observed in our simulation vs. 
configuration rank R. The inset shows a log-log plot of the number of unique configurations, n, as a 
function of the number of realizations r. The straight line asymptote is the power law n = 2, 244r'^'°^^''. 
In the four most probable configurations (large dot at Rank = 1), which all have substantially the same 
probability, both axes /3 and 7 form Mobius strips (|t^| = \T,y\ == 2), the streamlines of axis a form a spiral 
(Aq = 0), and those of both axis (3 and axis 7 form saddle points {Ia,f3 = — 1, = A-y = 4). 



Listed below are the probabilities of occurrence in our database of a number of configurations 
that may be of special interest. These probabilities are not densities because they do not include 
the appropriate Jacobian. Berry and Dennis have presented a Jacobian appropriate for some 3D 
L line properties [12], but this Jacobian does not appear to be suitable for Mobius strips, r-rings, 
and a-cones in T,q. The reason is that these structures and their statistics can change importantly 
when the plane of observation S is tilted away from Sq by a small but finite, amount, a property 
not shared by the Jacobian in [12]. Because the probabilities listed below are not densities they are 
not easily measured experimentally; they are, however, amenable to calculation. Fig. 12 presents 
some examples of transformations of Mobius strips, r-rings, and a-cones under tilt of the plane of 
observation; these transformed structures differ importantly from those appearing in Sq. A more 
complete discussion of transformations under tilt will be presented elsewhere. 
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FIG. 12: Transformations of Mobius strips, r-rings, and a-cones. As the plane of observation S is tilted 
away from Sq, different sets of ellipses occupy the surrounding circle. The structures, indices, and statistics 
of the Mobius strips, r-rings, and a-cones generated by these ellipses can change importantly with tilt angle. 
For reasons given below, in the examples shown here the radius r of the surrounding circle is r — 0.1, instead 
of the value used throughout this report, r = 10~*; in all the cases shown, when the tilt angle is zero, the 
Mobius strips, r-rings, and a-cones are, except for a scale factor, virtually identical for r ~ 0.1 and r = 10~^. 
(a) The Mobius strip in Fig. [T|3 for zero tilt angle, i.e. E = Sq. As in Fig. [T|d, axis endpoints are shown by 
small black (white) circles if the endpoint lies above (below) the circle of ellipse centers (small gray circles). 
As can be seen, just like in Fig. [Tjs the endpoints and the ellipse centers form interlocking rings with four 
crossings, and r = +2. (a') The strip in (a) for E tilted away from Sq by 5°. Here the rings still interlock, 
but with only two crossings instead of four, and r — +1 instead of r = +2. (a") The strip in (a) for S 
tilted away from Sq by 10°. Here the rings no longer interlock, and t — Q. (b) The opened r-ring in Fig. 
[I]!' for a tilt angle of 5°. Here the axes oscillate through only one sinusoidal cycle instead of through two, 
and dr = -1-1 instead of +2 as in Fig. [5| (c) The a-cone in Figs. [l|; and [6] for a 5° tilt angle. Here the 
endpoint spiral has only one turn instead of two, and dr = +1 instead of -1-2 as in Fig. [6]:. The reason for 
using r — 0.1 instead of r — 10^^ is that the tilt angle required for a given transformation decreases with r, 
becoming so small for r = 10~^ that experimentation becomes impractical. For rg = 0.1, however, the tilt 
angle is several degrees, and we use the larger value of r to demonstrate that the transformations described 
here can be accessible to experiment. 



From our simulation we find the following: 

(i) As expected, positive and negative values for all indices appear with equal probabilities. 

{ii) The fraction of L lines or L line segments surrounded by Mobius strips is 0.359, the 
remaining lines are surrounded by r-rings. This fraction holds for both axis (3 and axis "f strips. 

{in) The fraction of L lines or L line segments in which both axis (3 and axis 7 of the surrounding 
ellipses generate Mobius strips (r-rings) is 0.185 (0.468). 

(i?;) The fraction of Mobius strips whose vro endpoint curves form simple closed curves is 0.815, 
the remaining endpoint curves are figures of eight. This fraction holds for both axis /3 and axis 7. 

{v) The fraction of r-rings whose tq or vro endpoint curves form simple closed curves is 0.718, 
the remaining endpoint curves form figures of eight. This fraction holds for both axes j3 and 7. 

(m) The fraction of axis /3 Mobius strips (r-rings) with Ip^^ = -|-1 and A^ = is 0.808 (0.547), 
for the remaining strips (rings) with /^^^ = -|-1, A^ = 4. The same fractions with subscripts (3 and 
7 interchanged hold for axis 7 Mobius strips (r-rings). 
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(vii) The fraction of Mobius strips (r-rings) with = is 0.408 (0.228), for the remaining 
strips (rings), Aq = 4. These fractions hold for both axis (3 and axis ■j Mobius strips (r-rings). 

(vii) The fraction of a-cones with = +1 and A^ = is 0.586, for the remaining = +1 
cones Aq = 4. 

Although there are an enormous number of other statistical combinations that are easily ob- 
tained from our database, we believe that the above are likely to be of the greatest interest since 
they relate to the major geometric properties of, and interconnections between, the different struc- 
tures. 

V. SUMMARY 

The axes of the ellipses surrounding lines of linear polarization, L lines, in a 3D optical ellipse 
field have been shown to generate three distinctly different types of structures. In a plane normal to 
the direction of polarization of an L point on an L line the major and minor axes of the surrounding 
ellipses that lie on small circles surrounding the point generate Mobius strips or r-rings (rippled 
rings), whereas the major axes of these ellipses generate cones (a-cones). The Mobius strips have 
two whole twists and can be either right- or left-handed screws. These three different types of 
structures were characterized by 15 indices: 12 winding numbers that characterize the projections 
of the different structures onto three orthogonal planes, Sq, tq, and ttq, and three indices that 
characterize the streamlines formed by the projections onto Sq. Analytical (numerical) results 
were given for all indices in Sq (in tq and in ttq). Binary and ternary selection rules were presented 
that reduce the number of possible configurations from 839, 808 to 8, 248, of which 5, 562 have been 
observed in a computer simulation; probabilities of occurrence were given for the most important 
of these configurations. 

Experimental measurement of the various structures require determination of all three orthogo- 
nal components of the optical field; such measurements are currently feasible both in the microwave 
[9, 10, 18, 19], and in the optical [20 — 27], regions of the spectrum. 
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